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We have studied the zero-temperature statistics of the charge transfer between the two edges of 
Quantum Hall liquids of, in general, different filling factors, i/q^i = l/(2mo,i + 1), with mo > mi > 0, 
forming Mach-Zehnder interferometer. General expression for the cumulant generating function in 
the large-time limit is obtained for symmetric interferometer with equal propagation times along the 
two edges between the contacts and constant bias voltage. The low- voltage limit of the generating 
function can be interpreted in terms of the regular Poisson process of electron tunneling, while 
its leading large-voltage asymptotics is proven to coincide with the solution of kinetic equation 
describing quasiparticle transitions between the m states of the interferometer with different effective 
flux through it, where m = 1-1- mo -I- mi. For m > 1, this dynamics reflects both the fractional charge 
e/m and the fractional statistical angle 7r/m of the tunneling quasiparticles. Explicit expressions for 
the second (shot noise) and third cumulants are obtained, and their voltage dependence is analyzed. 

PACS numbers: 73.43.Jn, 71.10.Pm, 73.23.Ad 



I. INTRODUCTION 

Electronic Mach-Zchndcr interferometer (MZI) 
can be realized with the edge states of the Quantum Hall 
liquids (QHLs). Together with the quantum antidots 
, MZIs in the regime of the Fractional Quantum Hall 
effect (FQHE) are expected W§\ to be useful for obser- 
vation of the fractional statistics of FQHE quasiparticles. 
In contrast to fractional quasiparticle charge, which has 
been confirmed in several experiments [13, there 

is no commonly accepted observation of anyonic statis- 
tics of the quasiparticles, which remains a challenging 
experimental problem. Currently, this problem attracts 
interest in the context of solid-state quantum computa- 
tion, since individual manipulation of anyonic quasiparti- 
cles involving their braiding provides an interesting possi- 
ble basis for implementation of the quantum information 
processing |lll4l3l |. However, in typical interferometer- 
based experimental set-up the quasiparticles emerge as 
a continuum of gapless edge excitations that should be 
described by a ID field theory [ij]- Individual quasipar- 
ticles can be realized in this theory only asymptotically 
in a special limit. In the fractional edge-states MZI, such 
a limit occurs at large voltages, when the system is char- 
acterized by the Hamiltonian dual to the Hamiltonian of 
the initial electron tunneling model of the MZI. The lat- 
ter is perturbative in electron tunneling at low voltages, 
and is much better defined, since weak electron tunneling 
is probably the most basic process in solid-state physics. 

In the dual description derived from the electronic 
model by the instanton technique [1, [l^ , MZI acquires 
m = niQ + mi + l different quantum states which differ by 
the the effective flux $ through it. This flux contains, in 
addition to the external flux ^exi a statistical contribu- 
tion earlier also found \W\ for the antidots. The tunneling 



of each quasiparticle changes this effective flux by ±$0 
and therefore switches the MZI from one flux state into 
another. Since quasiparticle also carries the charge e/m, 
the change of flux by ±$0 results in the change of the in- 
terference phase for the quasiparticles by 27r/m and the 
corresponding change of the rate of coherent tunneling 
through the interferometer. Summation over the m flux 
states in calculation of the physical quantities restores 
$0 periodicity of their $ea; dependence. At low voltages, 
this periodicity is guaranteed by $0 periodicity of the 
electron tunneling amplitude. 

In this work, we follow the approach to the MZI that 
allows us to obtain the uniform description of its trans- 
port properties in both regimes of electron and quasipar- 
ticle tunneling. We consider the standard MZI geometry 
(Fig. 1) with two tunneling contacts between the two 
effectively parallel edges of QHLs, but allow for, in gen- 
eral, different filling factors, z/0,1 = l/(2mo,i -I- 1), with 
iTiQ > mi > 0, of these edges. In the symmetric case 
of equal propagation times along the two edges between 
the contacts, the corresponding ID field theory permits 
an exact Bethe ansatz solution [8|. Making use of this 
solution, we calculate the zero-temperature full counting 
statistics [13] of the charge transferred between the two 
edges forming MZI. The transferred charge distribution 
is shown to refiect the anyonic braiding statistics of the 
tunneling quasiparticles, and decomposition of electrons 
into quasiparticles with increase of voltage. 

The paper is organized as follows. In Section 2 we 
introduce the model of symmetric MZI and its Bethe 
ansatz solution. In Section 3, we use this solution to 
derive the general expression for the generating function 
P(^) of the distribution of the transferred charge. The 
method employed in this calculation generalizes to the 
MZI transport through two point contacts the method 
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FIG. 1: Mach-Zehnder interferometer considered in this work: 
two point contacts with tunnehng amphtudes Uj formed at 
points Xj, j = 1,2 between two co-propagating edges of QHLs 
with different filling factors fo and vi. The edges are assumed 
to support one bosonic mode each, with arrows indicating 
direction of propagation of these modes. 



developed earlier by Saleur and Weiss [l8| for a sin- 
gle point contact. The charge transfer statistics found 
for one contact [l^ demonstrated fractionalization of 
electron charges with increasing bias voltage across the 
contact. In the case of MZI, one can expect that in 
addition to charge fractionalization, braiding properties 
of anyonic quasiparticles should emerge in the charge 
transfer statistics. To see this, we analyze our result 
for the transfer statistics in different regimes. We find 
that the logarithm of the generating function separates 
into two parts, which for very different absolute values 
of the electron tunneling amplitudes in the two con- 
tacts, U2/U1 ^ 1, become identical with the cumulant- 
generating functions of the two separate contacts. In gen- 
eral, however, each "single-contact" term accounts also 
for the interference between the contacts. At low volt- 
ages, these single-contact terms are combined in such 
a way that the charge transfer occurs only in electron 
units, with the low- voltage asymptotic describing the reg- 
ular Poisson distribution of tunneling electrons. In Sec- 
tion 4, we develop analogous qualitative interpretation 
of the large- voltage asymptotics of the Bethe ansatz re- 
sult for the generating function. The interpretation is 
based on the m-state model of the quasiparticle tunnel- 
ing described above. Following the method of Ref. [2^, 
the cumulant-generating function is calculated from the 
kinetic equation governing the quasiparticle transitions 
in the basis of the m flux states. Direct comparison of 
this generating function with the leading asymptotics of 
the Bethe ansatz result reveals their coincidence under 
a special choice of m quasiparticle interference phases 
which are found as = (k + (m — 1)tt -t- 27rZ)/m, where 
I = 0, m—l numbers the flux states of the interferom- 
eter. The Z-dependent, statistical, part of these phases 
agrees with the expected anyonic statistics of the quasi- 
particles, while the common phase is given by the elec- 
tron interference phase n with an additional phase shift 
ir/ni for even m. The equivalence of the two distributions 
in the large-voltage limit proves that the Bethe ansatz 
construction we implemented indeed describes the sta- 
tistical transmutation of the effective flux through the 



MZI. In Section 5, we use the obtained generating func- 
tion to find the first three cumulants of the transferred 
charge distribution. We study, in particular, the volt- 
age dependence of the first and the second cumulants 
proportional to the average current and the shot noise, 
respectively. The ratio of the two, the Fano factor, is an 
experimentally observable [1, [l^ characteristics of the 
electron-quasiparticles decomposition. Finally, we con- 
sider the third cumulant, which determines the asymme- 
try of the transferred charge distribution around average. 
In Section 6, we look at the special case m — 2 which al- 
lows to obtain these cumulants in terms of elementary 
functions for arbitrary voltages. The results of this work 
are summarized in the Conclusion. 



II. MODEL OF THE SYMMETRIC 
MACH-ZEHNDER INTERFEROMETER AND 
ITS BETHE-ANSATZ SOLUTION 

We start our discussion with the electronic model of 
MZI (Fig. 1) formed by two single-mode edges with filling 
factors z^; = l/(2mi-|-l), / = 0,1. Electron operator of 
the edge I is expressed using the standard bosonization 
approach as 

■<pi = (L>/27rw/)i/26e'''^'^'''*^/^+'''''^I. 

Here 4>i are the two chiral bosonic modes propagat- 
ing in the same direction (to the right in Fig. 1), 
which satisfy the usual equal-time commutation relations 

(x) , 0p (0)] = Z7rsgn(a;) ^jp. The Majorana fermions 
account for the mutual statistics of electrons in differ- 
ent edges, and D is a common energy cut-off of the edge 
modes. The Fermi momenta ki correspond to the average 
electron density in the edges, while the operators of the 
density fluctuations are: pi{x,t) = {y'ui/2n)dx(t>i{x,T). 

In the symmetric case of equal times of excitation prop- 
agation between the contacts along the two edges, the 
two combinations of the bosonic operators 4>i that enter 
the electron tunneling terms of the two contacts can be 
expressed as the values at points xi^2 of the same right- 
propagating chiral bosonic field: 



0_(x) 



(1) 



The Langrangian describing electron tunneling in the two 
contacts can then be written as: 



{DUj/n) cos[X(f>^{xj) + Kj 



(2) 



where Uj and Hj are the absolute values and the phases 
of the dimensionless tunneling amplitudes. The products 
of the Majorana fermions ^1^2 were omitted from the 
Lagrangian ([2]), since they cancel out in each perturbative 
order due to charge conservation. The phases Kj include 
contributions from the external magnetic fiux $ea; and 
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from the average electron numbers Nq^i on the two sides 
of the interferometer: 

27r[($ex/$o) + (A^o/i'o)-(A^iM)]+const = -k. 

The factor A = \/2m in the Lagrangian © foUows from 
the normahzation of the bosonic field which in the 
absence of tunneling is a free right-propagating chiral 
field. This field undergoes successive scattering at the 
two contacts by the tunneling terms of the Lagrangian. 
The scattering breaks the charge conservation and there- 
fore creates tunneling current. The applied voltage can 
be introduced at first as a shift of the incoming field of one 
of the edges: 0o — y/T^Vt. As one can see from Eqs. ((T|), 
such a shift translates into the shift of the tunneling field 
</)_ - Vt/\. 

The thermodynamic Bethe ansatz solution of the tun- 
neling model was developed 19] for a single-point tunnel- 
ing contact with — 2m by application of one-particle 
boundary S'-matrices (20| to a distribution of the bosonic 
field excitations (kinks, antikinks, and breathers) intro- 
duced through the massless limit of the "bulk" sine- 
Gordon model. This solution was generalized [8] to the 
two tunneling contacts relevant for the MZI problem by 
successive application of two boundary 5-matrices to the 
same distribution of the excitations of the model. For the 
charge transport, only the kink-antikink (and vice verse) 
transitions are important, and their boundary 5-matrices 
are written as 



{ak/T. 



— — _ 

-1 ^ I 



ai(ak—Kj) 



iiak/TjE)""-^ 



(3) 

Here the standard energy scales Tjb, j = 1,2, are used 
to characterize the tunneling strength at the individual 
contacts, and 

2v7r(i/[2(i - 1/)]) 

^W[2(l-^)]) ■ 

The explicit relation between the energy scales Tjb and 
electron tunneling amplitudes is given below. 



III. CUMULANT-GENERATING FUNCTION 
FOR THE CHARGE TRANSFER DISTRIBUTION 

At zero temperature, dynamics of the liquid should 
be described with only one type of quasiparticles, e.g., 
kinks, which fill out all available states with the "bulk" 
distribution p{k), with k being the quasiparticle momen- 
tum, up to some limiting momentum A defined by the 
applied voltage. Each momentum- fc quasiparticle under- 
goes successive scattering at the two tunneling contacts 
independently of other quasiparticles. The overall scat- 
tering process is described by the product of the two 
boundary S'-matrices given by Eq. ^ . Our goal is to find 
the cumulant-generating function InP(^) of the charge 
transfer between the two branches of the interferometer, 



which is defined, as usual, as the logarithm of the Fourier 
transform of the probability distribution function of the 
transferred charge. We measure the charge in units of the 
elementary electron charge by setting e = 1. The inde- 
pendence of the scattering events of quasiparticles with 
different momenta implies then that InP(^) can be found 
as a sum of logarithm of the generating functions of indi- 
vidual momentum states, and its long-time asymptotics 



lnP(0=i/ dkp{k)\np{k,C). (4) 
Jo 

The generating function of one state with momentum k 



p{k,0 = l + rc{k){e'^ -I) 



(5) 



is defined by the total transition probability Tc{k) of 
the momentum-/c kink into antikink. Taking the prod- 
uct of the scattering matrices of the two contacts to 
find the total scattering matrix S2S1, and using the 
parametrization of the contact tunneling strengths as 
{TjB/cLf' = exp{dj/{m — 1)}, we can write the proba- 
bility Tc{k) from the corresponding matrix element of 
S2S1 as 

Tc(fc) = |(525i)-+p =B(T(^?2,fc)-r(^?i,fc)). (6) 

Here T{9j,k) are the transition probabilities in the indi- 
vidual point contacts: 

r(0,,fc) = \ST'+\^ = [l + /fc2(™-i)e-f.]-i, (7) 

and the factor B characterizes interference between the 
two contacts: 



B{TjB,n) = 



„2(m-l) _ „2(m-l) 
^2B ^IB 



(8) 



Without a loss of generality, we assume below a specific 
relation between the tunneling strength parameters of the 
two contacts, O2 > 9i, and write them as 61^2 — 0^ A^o, 
with A^o > 0. 

The aim of our subsequent derivation is to find the 
cumulant-generating function In P(^) in terms of the two 
generating functions In P5 for charge transfer in an in- 
dividual point contact that was found from the Bethe 
ansatz solution by Saleur and Weiss [l^. This deriva- 
tion does not need the explicit expressions for p{k) and 
A which can be found in [1£]. Following the approach 
for one contact, we first relate InP(^) in Eq. ^ to the 
effective tunneling current. To do this, we introduce the 
generalized tunneling probability Tc{u,k): 



Tciu,k) 



[l + {T^\k) 



l)e- 



(9) 



which is the solution of the following differential equation 
in the new parameter u: 



duTc{u, /c) = (1 - Tc{u, k))Tc{u, k) 



(10) 
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that satisfies the initial condition Tc{u,k)\u=o = Tc{k). 
One can extend Eq. ^ for p to include the parameter 
u through the substitution Tc{k) — ^ Tc{u,k). Equation 
© shows then that the logarithm of p extended this way 
can be expressed as 



lnp==ln[l + Tc(fc)(e"+^ 



1): 



I -S— 



(11) 



Calculating the derivatives of Eq. (fTTj) with respect to u 
and ^, and using Eqs. ([9]), one can see that 

— id^ \np = du Inp + tc{u, k) ~ tc{u + i^, k) . (12) 

Combining Eq. (fT2|) with (|4]) in which p is extended to 
include the parameter u, one sees that the cumulant- 
generating function satisfies the following relation 



a,4lnP(w,e)A= / dkp{k)Tc{u + iS,,k) = I{u + iS,,V), 
Jo 

(13) 

which expresses it through the total tunneling current 
I{u, V) in the two contacts that is defined by the gener- 
alized tunneling probability tc{u, k). 

As the next step, one substitutes Eq. ([6]) into ([9|) and 
casts the total tunneling probability tc{u, k) into the fol- 
lowing form 



Tc(u, k) 



^ e"sinhAgo ^^^ 
sinhA6l(M) ^ '^^ 



'±A6'(u),fc), (14) 



where A9{u) is defined by the conditions that 

cosh A9{u) = cosh A9a + S(e" - 1) sinh A9a , (15) 

and A9{u) > 0. Differentiation of Eq. (fT5|) shows that the 
coefficient in Eq. (jl4p in front of the sum can be written 
as the derivative of A9{u): 



duA9{u) = B 



e" sinhA6lo 
sinh A6'(w) 



(16) 



Using Eqs. (IT?)) and ([T5| in the definition of the tunneling 
current I(u, V) (jl3p we obtain an important relation 



liu, V) = duA9{u) J2 ±h/n 
± 



'tA0(m),V^), (17) 



which expresses the derivative of the cumulant- 
generating function (1131) for the charge transfer in a sym- 
metric interferometer in terms of the tunneling current 
Ii/m in one point contact between the two edges with 
the effective filling factor v = 1/m. The current in one 
contact has been calculated [l^ from the Bethe ansatz 
solution, and its tunneling conductance 

Gi/„AV/Tb] = Gi/„[(l//a)e-^/[2(™-i)]] ^ i^^^^^o,V)/V 

is given at zero temperature by a universal scaling func- 
tion expressed in the form of the low- and high-voltage 
expansion series. Integrating Eq. (1171) over u, and using 



the result in Eq. (|13p , we express the generating function 
InP(^) ~ lnP(u,^)|„=o in the following form: 

.e-Ae(»5) rS+Aeac) 
lnP(0 = -Vti^ / + jd9 

•Gi/™[(We-''/''^""'^']- (18) 

The explicit expansion series for Gi/^ in Eq- ifTS]) allow 
integration in each order. The integration transforms 
the generating function (fT5)) into the sum of the two gen- 
erating functions InPg for charge transfer in individual 
contacts and gives our key result: 



lnP(0= ^^Ps{VlT,B.e 



(-l)^(A0o-Ae«))^ 



(19) 



The single-contact generating function In P5 is known in 
terms of the low- and high- voltage expansion series [ISj : 



lnPs(s,e'^) ^ c^i(m)^2n(i-i)(g™e _ 



anVt ^ — ' mn 

n—1 



s < e 



2n{i/— 1) — mi^^ 



1) , s> e 



A 



1 (n -t- 1)1 (3/2 + [1^ — l)n) 



where 



„A 



and (To is the conductance quantum. 

Equation (|19p representing the total generating func- 
tion In P as the sum of the two single-contact generating 
functions In P5 seems to suggest that the charge transfer 
in the MZI is divided into two independent processes as- 
sociated with the two point contacts of the MZI. Indeed, 
one manifestation of this division is that dependence of 
each of these processes on the bias voltage V (through 
the function P5) is determined, as in Eq. ^2^ . only by 
the characteristics energy scale T^b of the corresponding 
junction. Such a division, however, is not complete. The 
total generating function In P depends also on the charge 
dynamics in the interferometer as a whole, since each in- 
dividual charge transfer process in one contact triggers 
multiple charge transfers involving interference between 
both contacts of the interferometer. Information about 
this charge dynamics enters Eq. (fTO)) through the function 
AQ(u) determined by Eq. (fT5)) and sensitive to both in- 
terferometer contacts. Nevertheless, the interference can 
become irrelevant if the two contacts are strongly asym- 
metric. For T2B S> PiB, the generating function InP is 
well approximated at low voltages by the single-contact 
InPs defined by Pis, i.e. by the strongest electron tun- 
neling amplitude \J\ - see Eq. ([24]) below. At large volt- 
ages, the generating function In P is approximated by the 
single-contact In P5 defined by P2B, i-e., by the weakest 
electron tunneling amplitude U^- 
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A. Low-voltage behavior of the generating function Tjb are expressed through the electron tunneling ampli- 



At small voltages, when V < Tjbb'^ for both j = 1, 2, 
the low- voltage expansions of the single-contact generat- 
ing function In Ps (HOI) for both terms in Eq. ([T^ can be 
combined as follows: 

oo , s 



z — ^ mn ^ — ^ ^ 



mn 

n=l 3 = 1,2 



JB) 



r COsh{nAe{iO) . 
■ L cosh(nA0o) i - ^ ' 

Since cosh(nA0(i^)) is a polynomial of cosh(A6'(i^)), and 
the latter, according to Eq. (|15l) . is a linear function of 
e'^ , this expansion of the MZI generating function shows 
that at low voltages, the charge transfer between the 
edges of the MZI is quantized in units of electron charge. 

More explicitly, using the standard expansion of 
coshnx (see Eq. 1.331.4 in Ref. [2ll ) and the relations 
that follow from Eqs. © and P^ : 



cosh A6'(iO = cosh A6io[l + R{e'^ - 1)] 



\rj-im — l 

R=B tanhAOn = ' 



p. 



2(m-l) 
IB 



P. 



2{m-l) 



(22) 



2B 



cosh Af?o 



^2(m-l) „2(m-l) 
^IB +^2B 

2(riBr2B)(™-i) ^ 



we bring Eq. ([?T|) into the following form: 
c„(m) 



^nPio--ovtY^^^[j:{v/p,Bf-'r 



n=l 



mn 



V2] 



{[l + i?(z-l)r + nE 



(-i)'cr_7 



Z(2coshA6'o 



[l + i?(z-l)]"-''} 



21 
z=l 



(23) 



Equation ([25)1 quantifies our previous conclusions 
about properties of the MZI charge transfer statistics 
in the low- voltage regime. In the limit of strongly differ- 
ent contacts, P2B ^ Tib, one finds that R ^ 1 and 
coshA(?o ^ 1, so that the charge transfer statistics ([25)1 
approaches that of one point contact [l^ characterized 
by PiB (i.e., the strongest electron tunneling amplitude 
Ui) with corrections in P\b/T2b also quantized in the 
electron charge units. The nth term in the expansion of 
this statistics in powers of bias voltage V corresponds in 
this case to tunneling of exactly n electrons. By contrast, 
the nth order term of the general MZI transfer statistics 
(1231) involves transfer of all numbers of electrons up to n. 
In the lowest order in V, the MZI statistics reduces to 
the Poisson distribution, with the coefficient in front of 
(z — 1) in Eq. (|23p is equal to the average electron tun- 
neling current. One can check this starting from Eq. ^ 
by direct perturbative calculation , if the energy scales 



tudes Ui 



P,B = 2D[T{m)/Uj] 



l/(m-l) 



(24) 



B. Large-voltage behavior of the generating 
function 



At large voltages, V > PjBe'^, the combination 
large-voltage expansions of both single-contact generat- 
ing functions In Ps in Eq. (IT^ brings the total MZI gen- 
erating function into the following form: 

InP(e) = aoVt 2^1^1M ^ (l//T,-5)2"(i-")/- 

n=l j 

. ^ cosh(nA6'(^0/TO) _ . 
cosh(nA6lo/m) J ' 

In terms of the parameters introduced in Eq. (|22p. this 
equation can be rewritten as: 



lnP(0 = aoF.f:^^[E(?^.B/V^)^<-- 



E [H- R{z - 1) ± ([1 + R{z - 1) 



cosh"-" A6ln ) " 



(26) 



We can see again that in the asymmetric limit, P2B ^ 
Tib, when R ^ 1 and coshA^o > 1, Eq. ^ for the 
charge transfer statistics reduces to that of a single point 
contact (isj . The dominant contact is now character- 
ized by the larger quasiparticle tunneling amplitude W2 
(i.e., smaller electron tunneling amplitude U2) and cor- 
responding energy scale r2B, related to W2 as: 



PjB = 2mD[Wj/r{l/T 



(27) 



The ?ith order term in the expansion ((26|) of the gen- 
erating function corresponds to the transfer of the frac- 
tional charge n/m by n quasiparticles. One can see, how- 
ever, that the MZI transfer statistics ([26)) does not con- 
tain in general the terms e'^/™ that would correspond 
directly to transfer of individual quasiparticles of charge 
1/m. In particular, the n = 1 term of the expansion 
that gives the leading large-voltage contribution to the 
statistics, can not be interpreted as a Poisson process of 
tunneling of independent quasiparticles, in contrast to 
the leading low-voltage term that did represent Poisson 
process of individual electron tunneling events. The rea- 
son for this is the m-state dynamics of the effective flux 
through the interferometer associated with the quasipar- 
ticle tunneling, which introduces correlations in the tun- 
neling process. These correlations can be most easily un- 
derstood in the description of the quasiparticle tunneling 
based on kinetic equation. Such an equation is discussed 
in the next Section. 
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IV. KINETIC EQUATION FOR THE 
LARGE- VOLTAGE CHARGE TRANSFER 

To derive kinetic equation that reproduces the large- 
voltage asymptotics of the generating function (pS)) . and 
therefore provides a simple physical picture of the dy- 
namics of quasiparticle tunneling, we start by rewriting 
this asymptotics in terms of the quasiparticle tunneling 
amplitudes Wj. From the last relation in Eq. ([2^ . we 
have 



cosh(A0o/"^) = - 



'T2B' 



2B - 



(28) 



Using this equation to transform the leading, n — 1, term 
in Eq. (|25|) . and replacing the energy scales TjB by the 
quasiparticle amplitudes Wj with the help of Eq. (P?]) . 
we have: 



InP(z) = tK{V) [2W1W2 cosh 



where X(F) = cro^(2mL>/F)2(™-i)/"ci(l/m)/r2(l/TO). 

Kinetic equation describing the quasiparticle tunneling 
can be written down based on the following considera- 
tions. The general picture of the quasiparticle dynamics 
in the MZI discussed in the Introduction implies that the 
quasiparticles create statistical contribution to the effec- 
tive flux through the interferometer. Because of this sta- 
tistical contribution, the MZI can be found in m separate 
states which differ by the effective flux, with each tun- 
neling quasiparticle changing successively the state I into 
I — 1 modulo m. Since the total rates of the quasiparticle 
tunneling in the MZI depend on the effective flux, such 
dynamics of flux introduces correlations into quasiparti- 
cle transitions, separating naturally the processes of suc- 
cessive quasiparticle tunneling events into the groups of 
m transitions. As usual, to cast the kinetic equation gov- 
erning this flux dynamics into the form appropriate for 
the calculation of the cumulant-generating function, we 
multiply transition probabilities by a factor ^i/™ = e*?/™ 
that keeps track of the transferred charge. Then, we in- 
troduce the probabilities d;,„(i) that at time t the MZI 
is in the state I and n quasiparticles have been trans- 
ferred through it. Combining these probabilities into an 
TO-dimensional vector Qi{z, t) = ^„ di^n{t)z"-^™ , one can 
write the kinetic equation in the following matrix form: 



dtQi{z,t)^J2^'^(^)u'Qi'iz.t)- 



(30) 



According to the qualitative picture of quasiparticle 
tunneling discussed above, the transition matrix has a 
simple form, with the only non-vanishing elements are 
those on the main diagonal, 1 = 1', and those with I = 



M{z) 



Ll' 



-ll 



(31) 



Here the Kronecker symbol 5i^ii is defined modulo m. The 
leading large-time asymtotics of the generating function 



for the probability distribution evol ving according to the 
kinetic equation ([5(7| is (see, e.g., [23]): InP(z) = th., 
where A is the maximum eigenvalue of the transition ma- 
trix (PT|) . The structure of this matrix shows directly that 
the characteristic equation det(M — A) = has the form: 



[] (7( + A) - z Jl 7i 

1=0 



(32) 



i=0 



The maximum eigenvalue A is the solution of this equa- 
tion which goes to zero when z — > 1 , since all other eigen- 
values of the matrix M{z = 1) are negative. 

Before trying to establish the general relation between 
the generating function obtained from this equation and 
the generating function , we consider the simple case 
m = 2. In this case, Eq. ([29]) for the large- voltage asymp- 
totics of the generating function (l^5t . can be simplified 
further. First, we have from Eq. (l22t : 



(29) cosh( 



Ae'(z) 



R, 



2 ) = [cosh^(^) + -(z-l)coshA0o 



1/2 



Using this relation, Eq. (gl]), and Eq. (US]) with m = 2, 
we transform Eq. ((29)) into: 



InP(z) = tK{V)[ W]f + \Wl + Wle'^p 
3 



On the other hand, in the kinetic-equation approach, 
Eq. (j32p be solved readily for m = 2 giving the following 
expression for A: 

A = ( [( E ^^0 ' + 47071 (^ - 1)] - E 7.) /2 • (34) 



This equation describes the generating function for the 
statistics of any transfer process consisting of the two 
steps with the rates 70,1, e.g. , incoherent charge transfer 
through a resonant level [22] ■ Comparison of Eqs. ([M)) 
and ([55]) shows that the generating function ob- 
tained from the kinetic equation reproduces the large- 
voltage asymptotics ([55]) of the generating function (|25p , 
if the tunneling rates are taken as 



70 



1= K{V)\Wi + W2e"*''-'\^, 0; = (K + 7r)/2 + 7rZ. (35) 



Equation (I35p for the tunneling rates agrees with the 
physical picture of quasiparticle tunneling discussed 
above. Statistical contribution to the effective fiux 
though the MZI means that tunneling of each quasipar- 
ticle changes the phase between the interferometer con- 
tacts by 27r/m = tt in agreement with the quasiparticle 
anyonic exchange statistics. Equation (1351) also shows 
that the quasiparticles see a phase shift 7r/2 in addition 
to the phase k/2 induced by the external magnetic flux. 
The origin of this phase shift is discussed below. 
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Following this logic, we look for solution of Eq. (|32|) in 

the case of arbitrary m taking the tunneling rates 7; as: 

-fi= K{V)\Wi+W2e'''>"\'^, ^ (t) + 27rl/m, (36) 

with some unknown (/). With these tunneling rates, 
Eq. dsn) reads 



_[J_(cos0H ^ ^ 

1=0 



2W1W2 



n 

1=0 



11 



2K{V)WiW2 



(37) 

where A = A/K{V). To further transform Eq. we 
use the basic identity: 



X 



X — e 



227r/ /rn \ 



(38) 



1=0 



and the two identities that follow directly from it: 

n'"^ 0/ r, , m-l 
9. ens — = 2 COS( — m + TT ^ ) , 



m — 1 



2™-i TT [cos(/); + cosh( — )] = cosh Aff (39) 



(-)'"cos(m(/ 



The first one is obtained essentially by taking x = — e*"^ 
in Eq. ([55]) . while the second one follows from the first if 
the sum of the cosines is transformed into their product. 
By direct comparison of the second identity in ((39|) with 
Eq. dSll) we see that A as defined by lnP(z)/i in Eq. ^ 
indeed solves Eq. if 



coshA6» = (-)'"cos(to0) + I T1 



11 



2 K{V)WiW2 



(40) 



Making use of Eq. (|38)) one more time, we calculate the 
product on the right-hand-side of Eq. (^0)) : 



m — 1 

n 

/=0 



11 



K{V)WiW2 



(41) 



One can see that with this expression for the prod- 
uct, Eq. (|40| precisely coincides with the definition of 
coshA0(z) by Eq. ([25]). Indeed, combining all three re- 
lations in Eq. p2|) one can express coshA0(z) as 



coshA^(z) = -C0SK+ ' If ^f, ^z. (42) 

Replacing T^b in Eq. (H^ with the amplitudes 
through Eq. (P7)) . we see that Eq. precisely coincides 
with Eqs. (jlo]) and (gl]), if cos k = (-)(™+i) cos{m(j)), i.e., 
if the phase </> is chosen to satisfy the condition 



m(j) = K + (m — l)7r . 



(43) 



The deviation of the phase (j) from k/tti in this equation 
is important only for even m. In this case, it produces 
the shift (m — l)TT/m of the interference phase in the 
quasiparticle tunneling rates from the value induced by 
the external magnetic field. This shift coincides with the 
phase acquired by one of the two Klein factors of the MZI 
quasiparticle tunneling action of the MZI [T^ in the fiux- 
diagonal representation due to the m-power periodicity 
condition for the Klein factors, which corresponds physi- 
cally to the requirement of the proper exchange statistics 
between electrons and quasiparticles. This phase shift en- 
sures, actually, that there is no shift in the interference 
pattern of the tunnel current in the interferometer (see 
the final results below) between the regimes of electron 
and quasiparticle tunneling. 



V. CUMULANTS OF THE CHARGE 
TRANSFER DISTRIBUTION 

So far, we have established the interpretation of the 
low- and high- voltage asymptotic behavior of the charge 
transfer statistics in terms of, respectively, tunneling of 
individual electrons and quasiparticles. In this Section, 
we calculate the charge transfer cumulants in these two 
limits, with the emphasis on the quasiparticle limit which 
exhibits the non-trivial behavior of the cumulants. We 
also will use the generating function found in Sec. IIIII 
for arbitrary voltages to calculate the full voltage de- 
pendence of the cumulants, and to study the crossover 
between the two asymptotic regimes of electron and 
quasiparticle tunneling. The cumulants of the charge 
N{t) transferred through the interferometer during a 
large time interval t can be found from the cumulant- 
generating function (jl3p by the standard relation: 

{N^{t))Jt - di lnP|„=oA = 5r i/(u, V)\u=o , (44) 

where the current I{u, V) is given by Eq. pT|) . 

The first cumulant gives the average tunneling current: 
I{V) = I{0,V) = {N{t))/t. At arbitrary vohages, the 
average current was calculated before in [15]. Its large-T^ 
asymptotics can also be obtained directly from Eqs. ([55)) 
and ([57)) : 



I{V) 



1 



(45) 



where the interference factor B ([S|) can be expressed in 
term of the quasiparticle amplitudes Wi,2 as 



B{W,,k) 



-I- VFa^e 



m^iK 1 2 



(46) 



The second equality in Eq. (j45p provides direct interpre- 
tation of the asymptotics of the average current in terms 
of the quasiparticles transitions Q • It can be proven for- 
mally by making use of the identity 



Wi 



W2d- 



2UW2 



2K{V){Wl - Wi) 



(47) 
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which can be obtained by directly differentiating indi- 
vidual rates 7; in this equation. On the other hand, 
differentiating the product of all "fis together, as given 
by Eq. (HJ), and using Eq. (|43p . we obtained the sec- 
ond equality in (j45p . This equality agrees naturally with 
the simple solution of the quasiparticle kinetic equation, 
which gives the average tunneling rate as inverse of the 
average tunneling times in the different states of the in- 
terferometer. 

The second cumulant (iV^(t)) defines the spectral den- 
sity of the current fluctuations at zero frequency 5/(0) — 
{N'^{t))c/t, which at zero temperature reflects the shot 
noise associated with the charge transfer processes. From 
Eq. (|17|) , one can write the first derivative of the current 
as 



duI{u,V) = I{u,V) X du\n[duA0{u)] 



(48) 



Substituting this formula into Eq. (pi)) and calculating 
the derivatives of A9{u) for w — we obtain ex- 
pression for the spectral density of current: 



5/(0) (l-ScothAe'o)/-B- 



,ie,,V). (49) 



It is convenient to characterize the short noise repre- 
sented by this spectral density through the Fano factor 
F defined as = ^/(O)//. In the case of MZI, the Fano 
factor reflects both the charge and statistics of the tun- 
neling excitations and illustrates the transition between 
the electron and quasiparticle regimes. In the low-voltage 
limit, = 1 as a result of the regular Poisson process of 
electron tunneling. To find the Fano factor in the quasi- 
particle, large- voltage, limit, we start with Eq. (|49l) which 
gives the following general expression for F: 

F = 1-b{ coth A0O - J2 deh/Ud^^V) 

j 

•[^(-l)^/i/„,(0„F)]-'}. (50) 
j 

Using Eq. (1^^ . the fact that in the large- voltage limit 
only one quasiparticle tunneling term oc W"^ can be kept 
in the current Ii/m-, cf. Eq. (j45p . and that with the 
parametrization of the energy scales Tjb with 9 intro- 
duced above, ^^^"^ 



W'^ (X e^/™, we get from Eq. ([50| 



F ^l-B' 



2m 



1 



1^2 



}• (51) 



The Fano factor ([5T|) corresponds to the dynamics of 
quasiparticle tunneling as described by the kinetic equa- 
tion ([50)1 . This can be seen by following the steps similar 
to that taken above for the average current. Applying the 
differential operator from Eq. (|47p to individual terms in 



the sum of the inverse tunneling rates 7;, one obtains 
directly the following identity: 



Wid 



Wi 



2<JW2 



2K{V){W( 



E 



7; 



On the other hand, replacing the sum of inverse 7;s in 
this equation with the corresponding expression from 
Eq. gll): 

Y,ir' = ^l[BK{V){Wi-Wl)], 



and performing differentiation, we see that the large- 
voltage asymptotics ([5T|) of the Fano factor can be writ- 
ten in terms of the tunneling rates 7; as 



^ = E7rV(E7r^)' 



This result agrees with the calculation 23] based directly 
on the kinetic equation. Because of the complex nature 
of the quasiparticle tunneling dynamics characterized by 
TO different tunneling rates 7/, -F is not equal simply to 
the quasiparticle charge l/m but varies as a function of 
parameters, e.g. the interference phase k, between 1/to 
and 1. 
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FIG. 2: The zero-temperature Fano factor F of the tunnel 
current in the Mach-Zehnder interferometer formed by two 

V = 1/Z edges, i.e. for m = 3, as a function of the bias voltage 

V for different degrees of asymmetry of the tunneling strength 
of the two contacts characterized by the T\b /T2S ratio. The 
solid curves corresponds to the case of complete constructive 
interference, n = Q; for the dashed curves, n = n. \n the latter 
case, _F = 1 identically for identical contacts, Tib/T2b = 
1. The curves illustrate the transition between the electron 
regime f = 1 at small voltages to the quasiparticle m-state 
tunneling dynamics at large voltages. The transition region 
is characterized by the Fano factor F reaching the minimum 
below the quasiparticle minimum l/m = 1/3. 

At arbitrary bias voltage V, the Fano factor F should 
be plotted numerically. Figure 2 shows F in the case 
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TO = 3 which corresponds, e.g., to tunneling between the 
two = 1/3 edges. The curves are shown for differ- 
ent degrees of asymmetry between the two contacts of 
the interferometer and two values of interference phase, 
maximum constructive interference, k = 0, and complete 
destructive interference, n — tt. The range of variation of 
F with the interference phase k decreases with increas- 
ing junction asymmetry. In general, the curves show 
the transition between electron tunneling with F ~ 1 
at small voltages V to quasiparticle tunneling at large 
voltages. In the quasiparticle regime, F is still can be 
significantly different from 1/3 because of the non-trivial 
three-state flux dynamics of the MZI. In particular, for 
identical junctions, the three total quasiparticle tunnel- 
ing rates p6p satisfy the relation: 7; c>c cos^ / 2) . Tak- 
ing into account Eq. (H^]) for the non-statistical contri- 
bution to the interference phase, one can see that under 
the condition of destructive interference, k ~ tt, the tun- 
neling rate in one of the flux states of the interferometer, 
Z = 0, is much smaller than the rates in the two other 
states. This mean that on the relevant large time scale 
set by the smallest rate, the three quasiparticles transi- 
tion that transfer interferometer from state I = back 
to itself happen almost simultaneously, so that the three 
quasiparticle charges 1/3 are effectively transferred to- 
gether, restoring F back to 1. 

Finally, we study the third charge transfer cumulant 
that characterizes the asymmetry around average of the 
transferred charge distribution, and has been measured 
experimentally for electron tunneling in metallic tunnel 
junctions - see, e.g., [llillll, and in quantum point con- 
tacts [2^. As for the other cumulants, the large-time 
asympotic of the third cumulant is linear in time, and it 
can be characterized by the coefficient C3 = {N^)c/t. To 
calculate this coefficient, we first find the second deriva- 
tive of the tunnel current from Eq. (|48]) : 

dfjiu, V) = (dl \nduA9{u) + [du ln9„A0(u)]2)/ - 

^ [3d^A9dlA9dg T {d^Aefdl]ly^{e T AO, V). (52) 
± 

Derivatives of A9{u) here can be found from Eq. (IT5|) . In 
particular, the coefficient in front of / in Eq. (|52p can be 
expressed as: 

dl\aduAe{u) + [du\nduA9{uf = 1 - 3cothA6'(u) 

xduA9{u) + [3 coth^ A9{u) ~ \\(duA9{u)f .(53) 

Substitution of Eqs. dM]), and (HH) into Eq. ^ 

gives us the coefficient C3 of the third cumulant: 

C3 [1 - 3BcothA6io + B^(3coth^ A6'o - 1)]/ 

[3^2(1 _BcothA^?o)ae/l/„,(^^,,F) (54) 



i=i,2 



^[--yB''dlhi,^{9,,v)\. 



The ratio F3 = C3// has also been suggested [27| as 
a possible alternative to the Fano factor to characterize 



the charge of the tunneling particles. Indeed, in a Poisson 
process, F3 is equal to the Fano factor multiplied by the 
tunneling charge. Therefore, in the case of MZI, F-j, (|54|) 
reduces to 1 in the low- voltage limit, as a result of the 
regular Poisson electron tunneling. In the quasiparticle, 
large- voltage, limit, repeating the calculation similar to 
that leading to Eq. (|50|) . we can relate both factors as 
follows: 



F3 = 3F-2 + B2<^2 



1 



121^1™ 



Ira 



1? {W^'^-Wl 
3 (W^l + W^) (VFl™ + W^"^) 



2m\2 



TO (W^l 



}. (55) 



VI. CHARGE TRANSFER STATISTICS FOR 



For general m, the results for the cumulants of the 
charge transfer statistics in the MZI discussed above can 
not be presented in a finite analytical form for arbitrary 
bias voltages. The situation is simpler for m = 2, when 
the kink-quasiparticles of the "bulk" sine-Gordon model 
that provide the basis for the Bethe-ansatz solution of 
the MZI transport are the regular fermions (though car- 
rying charge 1/2), and their distribution p{k) in Eq. (|4]) 
is the Fermi-Dirac step-function. In practice, the to = 2 
regime should take place in the MZI formed by the two 
edges of different Quantum Hall liquids, with filling fac- 
tors I' = 1/3 and v = 1. Individual tunneling contacts 
of this type have been realized experimentally [28]. Us- 
ing the Fermi-Dirac property of the distribution function 
p{k) for TO = 2, one can obtain the cumulant generating 
function InP(^) either by integration in Eq. ((4]), or equiv- 
alently, by direct substitution into Eq. (fT8|) of the known 
single-contact tunnel conductance G1/2, which can be ex- 
pressed simply as 6*1/2(5) — cro[l — arctan(2s)/(2s)]/2. 
The resulting generating function is: 

\\iP{C) = cTi^Vt ^ {yj{u) aicia.\)[l/yj{u)] 

u—i^ 



-(l/2)ln[l + y (^): 



n=0 



(56) 



where yj{u) = yj(0) exp{(-l)J (A6'(u) - A6'o)/2}, and 
y,{Q)^T,Bl{2V). 

This generating function can be combined with 
Eq. ((44|) to calculate the cumulants of the transferred 
change distribution for to = 2 using the same steps as 
in the previous section. Alternatively, one can substitute 
Eq. ([4]) into (l44t and perform the integration directly. 
Below we briefly discuss the first three cumulants ob- 
tained in this way. The first cumulant gives the average 
tunneling current in the MZI as fis}: 



I = aoB(T2l2-Tih), 



(57) 



where I, = arctan(y/2rj), and F^ = ^,5/4 = OWf/ir 
are the characteristic quasiparticle tunneling rates in sep- 
arate contacts. Using the fact that K{V) = a^D for 
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m = 2, one can see explicitly that the current agrees at 
large voltages with Eq. pst that follows from the quasi- 
particle kinetic equation with the tunneling rates (|35|) : 

_ ^do iri+r^e'-p _ 7071 

2 Ti+r^ "70 + 71- ^ ^ 

It is interesting to note that this agreement relies strongly 
on the shift of the interference phase (|43|) from the 
externally-induced phase At. 



rates l|35|) coincide, 71 =72, regardless of the relation 
between the individual rates Tj. Similarly to the case 
m = 3 illustrated in Fig. 2, the Fano factor reduces to 
electron value 1 even in the quasiparticle regime, if Fi ~ 
F2 and K — TT. In this case, one of the total quasipartical 
tunneling rates 7 is much smaller than the other, so on 
the relevant large time scale set by the smaller rate the 
quasiparticles effectively tunnel together, restoring F to 
1. 
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FIG. 3: The zero-temperature Fano factor F of the tunnel 
current for m = 2, i.e., in the Mach-Zehnder interferometer 
formed by edges with = 1/3 and v = 1, as given by Eq. (|59|1 . 
The solid curves corresponds to the case of complete construc- 
tive interference, k = 0; for the dashed curves, k = tt. In the 
latter case, F = \ identically for identical contacts. In gen- 
eral, the transition from the electron regime F = 1 at small 
voltages to the 2-state tunneling dynamics of quasiparticles at 
large voltages is characterized by the Fano factor F reaching 
the minimum in the crossover region. 

The second cumulant gives the following expression for 
the Fano factor at arbitrary voltages, including the tran- 
sition region between electron and quasiparticle tunnel- 
ing: 



F = 1 - 



iri 



■2FiF2(F2/i -ri/2 



Ti/i 
1 



r2/2 I (rf 



2(r?-rl)^.f^, ■'L^ + i + (y/2F,)2J 



}• (59) 



This equation is plotted in Fig. ^ and describes the 
transition from F = \ for electron tunneling at small 
voltages to 



1 



p i«|2 



2 (Fi+Fs) 



for quasiparticle tunneling at large voltages. One can see 
that the quasiparticle charge e/2 manifests itself most 
clearly for k = 0, when the total quasiparticle tunneling 
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FIG. 4: Alternative "Fano factor" F3 = Cz/I JBO]) related to 
the third cumulant of the tunnel current noise in the Mach- 
Zehnder interferometer with m = 2 and zero temperature, 
as a function of the interferometer phase k. The curves are 
plotted for several bias voltages V between the interferometer 
edges, and illustrate the transition between the electron and 
quasiparticles tunneling with increasing voltage. The tran- 
sition is characterized by the non-monotonous change of -F3, 
which reaches minimum at the intermediate voltages. 

Calculating the third cumulant, we find the following 
expression for the "alternative" Fano factor: 



F3 = - 2 



iri 



A{T,h-T^h){Tl-TlY ^^^^ 
3F, J, (Fj + lOF^F^, + 5F4 ) w 



F^ + 3r|, 



+ (-ir+i 



(y/2F, 



2(F?-ri) 
V 



(60) 
212 } ' 



where j' is defined by — 1,2, j' 7^ j. Equation 
([SO)) is plotted in Fig. |31 which shows F3 as a function 
of the interference phase n at several voltages. Voltage 
dependence of is qualitatively very similar to that of 
the Fano factor shown in Fig. ([3]): it approaches 1 at 
small voltages, in agreement with the underlying Poisson 
tunneling process of electrons. At large voltages, Eq. ([60)1 
reduces to the following form 



F.=ZF 



+ 4 (Fi 
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This expression can be understood in terms of the same 
two-state tunneling dynamics of the quasiparticles that 
was discussed above. Also similarly to the Fano factor, 
the voltage dependence of Ft, is non-monotonic, with a 
minimum between the regimes of electron and quasipar- 
ticlc tunneling. The main qualitative difference between 
the noise-related Fano factor and its third-cumulant al- 
ternative is that the minimum of can become nega- 
tive for some values of parameters (regime not shown in 
Fig. 15). 

VII. CONCLUSION 

Starting from the exact solution of the tunneling model 
of symmetric Mach-Zender interferometer in the FQHE 
regime, we have calculated the statistics of the charge 
transfer between interferometer edges. The obtained 
statistics shows the transition from electron tunneling at 
low voltages to tunneling of anyonic quasiparticles of the 



fractional charge e/m and statistical angle 7r/m at large 
voltages. Deep in the electron tunneling regime, the dy- 
namics of charge transfer is represented by the standard 
Poisson process. Dynamics of quasiparticle tunneling is 
more complicated and reflects the existence of m effec- 
tive flux states of the interferometer. The interference 
phase between the quasiparticle tunneling amplitudes in 
two contacts of the interferometer contains a contribution 
from the quasiparticle exchange statistics, making the 
quasiparticle tunneling rates in different interferometer 
states different. In general, the transition from electron 
to quasiparticle tunneling is reflected in the Fano factor 
F or its third-cumulant alternative F3 , which both reach 
minima in the transition region. However, in the regime 
close to complete destructing interference (interferome- 
ter phase K — n and equal tunneling strength in the two 
contacts), both F and F3 have electron value 1 for all 
voltages. 

V.P. acknowledges support of the ESF Science Pro- 
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